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We analyze the validity of the fluctuation-dissipation theorem for slow relaxation systems in the 
context of mesoscopic nonequilibrium thermodynamics. We demonstrate that the violation arises as 
a natural consequence of the elimination of fast variables in the description of a glassy system, and 
it is intrinsically related to the underlying activated nature of slow relaxation. In addition, we show 
that the concept of effective temperature, introduced to characterize the magnitude of the violation, 
is not robust since it is observable-dependent, can diverge, or even be negative. 



Many nonequilibrium systems in nature evolve in time 
following slow relaxation processes. Examples of this be- 
havior are usually encountered in glassy systems |jl| , poly- 
mers L granular flows || , foams [Q] , and crumpled ma- 
terials |5|] to mention just a few. A complete and satisfac- 
tory characterization of these systems constitutes nowa- 
days one of the most challenging issues of nonequilibrium 
statistical physics. The main feature of slow processes is 
that the relaxation time may exceed significantly the ob- 
servation time scale in such a way that the system can be 
considered as being permanently out of equilibrium. This 
peculiarity is the origin of a distinctive behavior which 
differs markedly from the case in which relaxation occurs 
in shorter time scales. The existence of aging regimes j(| 
and the violation of the fluctuation-dissipation theorem 
(FDT) constitute examples of this behavior. 

For all these reasons, the straightforward application of 
equilibrium concepts, appropriate to describe fast relax- 
ation processes, to out of equilibrium situations, inherent 
to slow relaxation dynamics, becomes in principle doubt- 
ful. However, our purpose in this Letter is to show that, 
when nonequilibrium thermodynamic concepts are ap- 
plied at the mesoscopic level H, one may justify many of 
the peculiarities of the behavior observed in glassy sys- 
tems. In particular, we will show that the violation of 
FDT is a natural consequence of the activated nature of 
the dynamics of a slow relaxing system. Starting from 
a more detailed description in which the system can be 
safely considered as near equilibrium and evolves via a 
diffusion process, we will show that the implicit elimina- 
tion of the fast variables, leads to an activated regime 
where the system becomes far from equilibrium and con- 
sequently the FDT is not fulfilled. Coarsening the level of 
description is then the origin of the violation of the FDT 
in strong glasses. Precisely, one way to characterize this 
violation is through the concept of effective temperature. 
We will also discuss the validity and robustness of this 
concept. 

It is well established that the evolution of many sys- 
tems can be described in terms of its energy landscape 



|g, [HJ, representing the (free) energy as a function of 
an order parameter or reaction coordinate 7 [jll). Com- 
plex systems exhibit a very intricate landscape with a 
great multiplicity of wells separated by barriers. Whereas 
at high temperatures the system may explore the whole 
landscape at low enough temperatures the dynamics re- 
duces basically to two elementary processes: a fast relax- 
ation toward the local minima via a diffusion process, and 
a slow activated process in which the system overcomes 
the barrier toward the next minimum. The presence of 
the barriers is thus the cause for the slow evolution of 
the system. Hence, the case of a single barrier captures 
the essential mechanism of the slow dynamics. To show 
how the activated nature of the slow evolution of the sys- 
tem can be responsible of some of the peculiarities of the 
response of glasses we will then focus on the simplified 
model of a bistable potential. 

It is then plausible to assume that the evolution of the 
system occurs via a diffusion process through its poten- 
tial landscape parameterized by the 7-coordinate, which 
will be characterized by the diffusion current J(7, t) and 
the corresponding chemical potential ^(7, t). As any 
diffusion process, it can be treated in the framework 
of nonequilibrium thermodynamics |j[2). Assuming lo- 
cal equilibrium in 7-space, variations of the entropy Ss 
related to changes in the probability density p{"f,t) are 
given through the Gibbs equation 



5s=-j; I f l( 1 ,t)S P (j,t)d 1 , (1) 



where T is the temperature. 

The entropy production inherent to the diffusion pro- 
cess, er = ds/dt, 
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follows from Eq. (Q), after using the continuity equation 
in 7-space, 
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(3) 



From that expression one then infers the relation be- 
tween current and thermodynamic "force" J (7, t) = 
— y J £(7, Y)-^jfJ-(Y : t)drj' . The assumptions of locality 
in 7-space, for which 1/(7,7') — -^(7)^(7 — 7')) an d 
ideality, imposing the form of the chemical potential 
fi (7, t) — ksTXnp (7, £) + $ (7), with fcs being the Boltz- 
mann constant, and $ (7) the bistable potential, provide 
the diffusion current in 7-space 

J (7, t) = _# e -*/^ JU^t (4) 

where D = ksL/p is the diffusion coefficient, taken to 
be constant as a first approximation. When this phe- 
nomenological relation is substituted into the continuity 
equation (|j) one obtains the diffusion equation 
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*(7) , (5) 



which governs the evolution of the average probability 
density. This result agrees with the one derived from a 
master equation [pH , and indicates that nonequilibrium 
thermodynamics can be used at mesoscopic level to pro- 
vide the fundamental kinetic laws of the Fokker-Planck 
type governing the dynamics. 

The probability density is subjected to fluctuations, 
which may be introduced through a random contribution 
to the current, J r (7,t), in Eq. (||) p3| , satisfying the 
fluctuation-dissipation theorem in 7-space 

(J r (7, t) J r (7', f )> = 2^ (P (7, *)> S^-^Sit-t'), (6) 

where (p(7,£)} is the solution of Eq. (||). The formula- 
tion of a FDT is intimately related to the fact that the 
system is in local equilibrium in 7-space. 

When the height of the energy barrier separating the 
two minima of the potential is large compared to thermal 
energy, which happens at low enough temperatures, a fast 
relaxation toward the local minima occurs, and the sys- 
tem achieves a state of quasi-equilibrium characterized by 
equilibration in each well. The chemical potential then 
becomes a piece-wise continuous function of the coordi- 
nates, m(7>*) = m(7i>*) (7o -7) + ^ (72,*) © (7 - To) , 
and consequently the probability density achieves the 
form 



p( 7 ,f) = p 1 (i)e- { * M ~*^ )}/fcBT e(7o-7) (7) 
+ p 2 {t)e-^-^l k - T Q{ 1 ~ lo ). 

Here p\ (t) = p(ji,t) and P2 (t) = p(j2,t) are the values 
of the probability density at the minima, (x) is the unit 



step function, and 7 , 71, and 72 are the coordinates of 
the maximum, and the minima of the potential, respec- 
tively. 

Hence, once the fast relaxation toward the local min- 
ima has occurred, the evolution of the system proceeds 
by jumps from one well to the other undergoing an acti- 
vated process. In this situation, a contracted description 
performed in terms of the populations at the wells can be 
adopted. This description corresponds to that of the two 
level model for a glass jl^, a minimal model which 
evolves according to an activated dynamics [|l6| confer- 
ring him the characteristic aging properties of glasses, 
closely related to hysteresis [[17| . Defining the integrated 
probability N(-f,t) — fl (x> d,'y' p{^' ,t), and by integration 
of the continuity equation (g) we obtain 



^-N(j,t) = -j s ( 7 ,t)-r( 7 ,t). (8) 

To proceed with the contraction of the dynamics from the 
diffusion regime to the two level regime we will introduce 
the integral operator P acting on a function in 7-space 
as 



P1K7) 



f 2 drfe*t^T 



d7e*/ fcBT V>(7)- (9) 



Projecting both sides of Eq. (j|) with P, using Eqs. (^) 
and (0), and evaluating the integrals using the steepest 
descent approximation, we obtain the equation governing 
the dynamics of the two state systemfDJ 
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(10) 



where ri\(t) = N(jo,i) and n 2 (i) = 1 — n\{t) are the 
"populations" at each side of the barrier. The value of 
the systematic current J(t) = PJ(j,t), which is the net 
current on top of the barrier, is given by 



J(t) = k—,ni — k^n-i = J_ 



(11) 



whereas J r {t) — PJ r (7, t), is the random current, whose 
correlation follows from Eq. (||) 

(J r (t) J r (t')) = (k^ (m) + k^ (n 2 )) 5 (t - t') . (12) 

In the previous expressions, k^ and £;<_ are the forward 
and backward rate constants 



k = iy$"(7i, 2 ) Hfr"(7o) 1 / $(7i,a) ~ *(to) 
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(13) 

It is important to highlight that Eq. ( p"2| ) evidences that 
the fluctuation-dissipation theorem is violated in the ac- 
tivated process. Only for fluctuations around equilibrium 
this equation becomes (J r (t)J r (t)) = 2k^n^ L 5{t — t 1 ), 
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which is the formulation of the fluctuation-dissipation 
theorem]!^]. In fact, Eq. (10) with this prescription con- 
stitutes a Orstein-Uhlenbeck process. The failure of the 
theorem, which was initially valid in 7-space, results from 
the coarsening of the description. When the dynamical 
description is carried out in terms of the reaction coor- 
dinate, the system progressively passes from one state to 
the other, which makes it possible the assumption of local 
equilibrium and the formulation of a mesoscopic nonequi- 
librium thermodynamics. However, when we describe the 
system in a characteristic time scale similar to the obser- 
vation time, we are only capturing the activated process, 
which is not near equilibrium and accordingly the FDT 
does not hold. 

The model we have introduced facilitates the analy- 
sis of the nonequilibrium response of the system. Let 
us consider, for example, the case of a dynamical ob- 
servable 0(7) (energy, density, magnetization, etc.). Its 
mean value, in the quasi-stationary regime, is (SO(t)) — 
J 0(j)5p{j,t)dj = (Oi - 2 )5n u where (SO(t)) = 
(0(t)) — (O) , and 0\ and O2 constitute the values of 
0(7) in the states 71 and 72, respectively. The response 
to an external perturbation —eO(j), plugged in at in- 
stant t w , will be characterized by the response function 
R{t,t w ) — d(SO(t)} /de(t w )\ e _> , This quantity can be 
calculated from Eq. (fl0|), yielding 
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where r = (fc_> + k<J) is the relaxation time for the 
activated process, which in view of Eqs. (13) is of the 
Arrhenius type. From Eqs. ( |ici| ) and (12) one can 
also calculate the correlation function jl9|], Co(t,t w ) = 
((0(t)-(0(t)))(0(t w )-(0(t w )))), which for t > t w 
and in the limit of large t, t w , is given by 



Co{t, t v 



(Ox 

-) (Snx(t w )) + k^n\ q } 



(15) 



R e J(t,t w ) = (Ox-Oz) 



eq 
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-{t-t w )/T 



(16) 



and is proportional to the time derivative of the equilib- 
rium correlation obtained from Eq. (|l5]), 



dt w 



(Oi - 2 ) 2 k^n e 1 q e 



-(t-t w )/i 



(17) 



We then recover the FDT relation Rq 



l/k B Td t „„Cf q 
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which holds irrespective the observable we are consider- 
ing. Out of equilibrium the FDT is not fulfilled, and 
its violation is usually quantified in terms of an effective 
temperature [Q, defined as 



Ro(t, t w ) 



k B T? tt dt, 



-Co(t, t v 



(18) 



-eff 



For the model we are considering, the effective tempera- 
ture, obtained from Eqs. dT3) and (flBl), becomes 



rpO 

1 eff 



1 



Ae-^/ T + (l - e-*™/ r ) 



(19) 



being A = *-<"iWKo>i-Oo) 

This expression reveals important conclusions. The ef- 
fective temperature does depend on the observable 
0(7) and explicitly on the waiting time t w . The depen- 
dence on the observable, which has also been found in 
a trap model for a glass J2(J and in experiments [pH , 
evidences that the effective temperature is not a robust 
quantity. Only for small deviation from equilibrium or 
when t w » r, one recovers the familiar result = T 
for all observables. It should be noted that our results, 
obtained by means of a non-mean field approach, differ 
from the ones following from mean field models |^, |22| 
(which yield an effective temperature independent of the 
observable) because the latter do not take into account 
the activated nature of the dynamics ||^, |J| . It also is 
worth to mention that, since T°jj depends on t w , the 
value of the effective temperature inferred from the slope 
of the FD plots, which represent the integrated response 
of the system x{t,t w ) — J. dt'R(t,t') against the cor- 
relation function, C(t,t w ), it is not the same as the one 
defined through Eq. (|lg|). 

Several interesting behaviors can be identified upon 
variation of the parameter A in Eq. ( |l9| ) . For < A < 1 , 
the effective temperature is higher than the temperature 
of the bath T, in agreement with the experimental mea- 
surements reported in ^7j. Contrarily, if A > 1 the effec- 
tive temperature is lower than the bath temperature T ; 
whereas, if A < 0, TQ* may diverge as predicted in |Q, 
numerically verified in |2j|, and experimentally suggested 
J2(J], or even become negative J2j|. All these cases are 
illustrated in Fig. [l], and arise from the peculiar behavior 
of the nonequilibrium response of an activated process. 
The effective temperature is essentially a measure of the 
ratio between the equilibrium and the nonequilibrium re- 
sponses of the system. When this ratio is smaller (larger) 
than one, then T° n < (>)T. A diver gence in Tf^ j occurs 
when the nonequilibrium response vanishes, and "nega- 
tive" effective temperatures would be caused by nonequi- 
librium responses having a different sign than its equi- 
librium counterpart. These anomalous behaviors can be 
tuned by a proper choice of initial conditions and observ- 
ables. To illustrate that fact, we have implemented our 
theory for two examples of bistable potentials: a quartic 
potential V(j) — 7 4 /4 + a7 3 /3 — 7 2 /2 — 127 , being a an 
adjustable parameter responsible for its asymmetry, and 



k^{n 2 (0))(O 2 -O ) 



At equilibrium, i.e. = J<_, the response reduces to in 
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FIG. 1: Qualitative behavior of T e ff/T as a function of t w /r 
for A < 0, < A < 1, and A > 1. 

the potential V(7) = — ecos(7) + 1/2 sin 2 (7), describing 
a monodomain magnetic particle |2Sj], Selecting differ- 
ent observables and initial conditions, we obtain in both 
cases the behaviors of the effective temperature shown 
in Fig. |]. In summary, we have shown that the origin 
of the violation of the FDT is the drastic elimination of 
variables one tacitly performs to model the system in the 
experimental time scale. At this level, the system evolves 
undergoing an activated dynamics, requiring big amounts 
of energy to surmount the barriers. Consequently, it is al- 
ways far from equilibrium and the FDT, a result strictly 
valid at or near equilibrium, is not fulfilled. In the more 
complete scenario, when instead of jumping between two 
states the system reaches a different state passing pro- 
gressively from intermediate configurations, i.e. diffusing 
in a configuration space, local equilibrium can be estab- 
lished. One can then proceed with the formulation of 
a mesoscopic nonequilibrium thermodynamics per- 
fectly compatible with the Fokker-Planck level of descrip- 
tion, whose underlying stochastic kinetics satisfies FDT. 
A widely-used way of quantifying the FDT violation is 
through the definition of an effective temperature. Our 
analysis shows that this concept suffers from a lack of ro- 
bustness, since its value depends on the dynamical vari- 
able we measure, and can diverge or even become nega- 
tive. All these problems limit the scope and question the 
usefulness of this quantity in the description of glassy 
systems where the activated dynamics is an unavoidable 
ingredient. 

The theory we have developed provides a useful frame- 
work to describe the behavior of systems with slow dy- 
namic bridging the macroscopic and the mesoscopic de- 
scriptions, by indicating the way to generalize local equi- 
librium concepts. 
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